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Bosonic Integer Quantum Hall effect in an interacting lattice model 
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We study a bosonic model with correlated hopping on a honeycomb lattice, and show that its 
ground state is a bosonic integer quantum Hall (BIQH) phase, a prominent example of a symmetry 
protected topological (SPT) phase. By using the infinite density matrix renormalization group 
method, we establish the existence of the BIQH phase by providing clear numerical evidence: (i) a 
quantized Hall conductance with \o-j;y\ = 2 (ii) two counter propagating gapless edge modes. Our 
simple model is an example of a novel class of systems that can stabilize SPT phases protected by 
a continuous symmetry on lattices and opens up new possibilities for the experimental realization 
of these exotic phases. 

PACS numbers: 73.43.-f, 03.65.Vf 


Introduction: Topological phases, which are quan¬ 
tum states of matter beyond Landau’s symmetry break¬ 
ing paradigm, are a main focus of modern condensed 
matter physics [1]. Well known examples are fractional 
quantum Hall (FQH) phases in two-dimensional electron 
systems [2] and quantum spin liquids that can arise in 
frustrated spin systems [1, 3, 4]. These examples fall into 
the category of intrinsic topological order [5], character¬ 
ized by their fractionalized quasiparticle excitations and 
the presence of long-range entanglement. More recently 
it was found that phases with only short range entangle¬ 
ment can also be topological non-trivial if certain symme¬ 
try is obeyed [6-9] . Such symmetry protected topological 
(SPT) phases have a gapped bulk with no fractionaliza- 
tion and are usually characterized by their anomalous 
gapless [10, 11] or intrinsic topologically ordered [12-15] 
edges. 

SPT phases in d-dimensional systems, protected by a 
global symmetry G, are generally classified by the coho¬ 
mology group [6, 7]. The first example dis¬ 

cussed in this context was the Haldane spin chain [16, 17], 
which is a one-dimensional SPT phase protected by either 
spin rotational (e.g., Z 2 x Z 2 ), time reversal or inversion 
symmetry [9]. The SPT phases protected by continuous 
symmetries can be described using effective topological 
field theories [12, 18-21]. For example, the bosonic in¬ 
teger quantum Hall (BIQH) state, which is a U{1) SPT 
phase (protected by the U{1) charge conservation), can 
be described by a mutual Chern-Simons theory [18, 19]. 
Along this line, BIQH states can be constructed using 
composite bosons in quantum Hall systems [19] and were 
found to be stabilized in two-component bosonic contin¬ 
uum models in a strong orbital magnetic field [22-24]. 

Bosonic SPT phases require interactions and have no 
simple non-interacting analogs, making it challenging to 
study and stabilize these exotic phases. Recently, several 
fine tuned exactly solvable lattice models have been in¬ 
troduced that realize SPT phases that are protected by 
discrete symmetries [10, 11, 14]. In contrast, the realiza¬ 
tion of more physical SPT phases, which are protected by 
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FIG. 1: (color online). The model is defined on a honeycomb 
lattice, and we assign time-reversal symmetry breaking back¬ 
ground flux in the system, (b) The unit cell is doubled to 4 
sites (along the y direction). The gauge held on the specihed 
link is taken as tt, tt — a/3. For the dashed links, the gauge 
helds are 0 (nearest neighbors), a/3 (next-nearest neighbors). 


continuous symmetries, has rarely been explored [25-27]. 
It is desirable to have concrete lattice models as they al¬ 
low for a relatively simple numerical study of topological 
phase transitions [28-30] and serve as a guide to where 
these phases might appear in nature, either in strongly 
correlated systems or optical lattices. Interestingly, an 
early study [31] (before the concept of SPT was intro¬ 
duced), has found a possible BIQH state in a Hofstadter 
type model by calculating the Hall coefficients axy = 2 
using exact diagonalization on small system size. How¬ 
ever, it is still unclear whether a C/(l) SPT phase is real¬ 
ized in such model, since the existence of anomalous gap¬ 
less edges, the most defining property of the SPT phases, 
has not been verified. 

In this letter, we introduce a simple lattice model for 
which we show numerically that it realizes a BIQH phase. 
Our model consists of hard-core bosons on a honeycomb 
lattice with correlated hopping which is subject to a back¬ 
ground gauge flux, and hence breaks time-reversal sym¬ 
metry explicitly. By using the infinite DMRG method 
[32, 33] on a cylinder, we obtain numerical evidence that 
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the model stabilizes two different BIQH phases with op¬ 
posite quantized Hall conductances a^y = ±2, as the 
background flux is tuned. Firstly, when adiabatically 
inserting 27r flux, we find that two bosons are pumped 
from the left edge of cylinder to the right edge, hence the 
ground state has a quantized Hall conductance axy = 2 

[34] , Secondly, using entanglement spectra as a probe 

[35] , we show that the ground state has two counter- 
propagating edge modes, and it fits theoretical expecta¬ 
tion perfectly. 

Model: Our model is defined on a honeycomb lat¬ 
tice with hard-core bosons at half filling described by the 
Hamiltonian 

H = - l)ajaj -k h.c. 

((*i» 

-)- Y, |^e*^'''(2n“ - l)blbi + h.c. 

{{ki)) 

(o> 

where a (6) is the annihilation operator for a hard-core 
boson on sublattice A (B) of the honeycomb lattice (see 
Fig. 1). The first two terms are the correlated hopping 
terms: the hopping of bosons on the sublattice A are cou¬ 
pled to the particle number of bosons on the sublattice 
B on the intermediate site, = bl.bk = 0,1; and vice 
versa. The last term is a standard nearest-neighbor hop¬ 
ping term. We also assign a background flux (j)^ (Fig. 
1(a)): the triangular plaquette with three vertices from 
the same sublattice has a flux tt — a, the small triangular 
plaquette with vertices from different sublattices has a 
flux a/3. Thus each hexagonal plaquette has a flux tt. 

When A = 0, the model has a 17(1) x 17(1) symmetry, 
which can be considered as either the particle conserva¬ 
tion on each sublattice, or the overall particle (charge) 
na + ny conservation plus the (pseudo) spin Ua — ny con¬ 
servation. However, the BIQH state only requires one 
17(1) symmetry to protect it. That is to say, we can al¬ 
low tunneling between sublattices (A ^ 0), so that the 
17(1) X 17(1) symmetry is broken down to a global 17(1) 
from overall charge conservation, without destroying the 
BIQH phase. We find that this is in fact true for our 
model: the BIQH state is robust as we turn on a finite 
A (< 0.7), and almost all the properties are the same as 
the A = 0 case. 

We use the infinite DMRG method [32, 33] to study 
the system wrapped around a cylinder with two different 
geometries (see supplementary materials) and find that 
the results do not depend on this choice. In the following, 
we set for simplicity A = 0 (unless specified otherwise). 
Due to the background flux tt (in each hexagon) we dou¬ 
ble the unit cell to 4 sites and we choose a gauge as shown 
in Fig. 1(b). Note that it is possible that certain choices 
of gauge and geometry will give rise to an extra flux for 
a non-contractible loop along the y direction (which can 



FIG. 2: (color online). Charge pumping in the BIQH state. 
Here we show two cases (the width is Wy = 8): A = 0.7, 
a = 7r/2, with Hall conductance a^y = 2; A = 0, a = —tt/S, 
with Hall conductance a^y = —2. For other system size, the 
results are almost the same. 

only be eliminated by a large gauge transformation). As 
the BIQH phase is gapped, such a global gauge flux will 
not cause any significant effect for our purpose. 

Based on the numerical simulations on systems of 
width Wy = 8 , 12,16 sites (corresponding to Ly = 2, 3, 4 
unit cells), we have obtained the phase diagram consist¬ 
ing two BIQH states with opposite Hall conductance axy- 
for a G (0,7r), axy = 2; for a G (— tt, 0), axy = —2. To es¬ 
tablish the existence of the BIQH phases, we now discuss 
in detail two of their characteristic fingerprints: (i) the 
ground state has a quantized Hall conductance \axy\ = 2 ; 
(ii) the ground state has two counter-propagating gapless 
edge modes. 

Quantized Hall conductance: The quantized Hall 
conductance is a hallmark of the quantum Hall state. In 
contrast to fermionic systems, the Hall conductance axy 
of a BIQH state is always quantized to an even number 
[19]. Numerically, we can use an adiabatic flux insertion 
to measure the Hall conductance axy- 27r flux insertion 
on a cylinder will pump axy particles from the left edge 
to the right edge of cylinder [34]. Such flux insertion 
can be implemented in the Hamiltonian by twisting the 
boundary condition in the infinite DMRG algorithm [36] : 
the bosons hopping around the cylinder pick up a flux 
The Hall conductance can then be written as [37]: 

p27T 

crxy= [d^JQ{%))]d<i>y, (2) 

^0 

(Q(<1>,)) = ^A($,)Q,($,), (3) 

i 

where Ai(<i)y) are the eigenvalues of the reduced density 
matrix when flux ^y is inserted, Qi{^y) is the corre¬ 
sponding ?7(1) quantum number. The numerical data in 
Fig. 2 clearly shows a quantized Hall conductance as 27r 
flux is inserted, and it is robust when turning on a finite 
A. 
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TABLE I: The degeneracies of energy levels of the entangle¬ 
ment Hamiltonian (edge mode) in two different charge sectors 
ANa + ANb = 0 and ANa + ANb = 1. 


Levels 

mode 

ky 

0-sector 

1-sector 

Ground state 

- 

0 

1 

2 

1st excited state 

charge 

271 / Ly 

1 

2 


Spin 

— 2'K/Ly 

3 

2 


charge 

4-KlLy 

2 

4 

2nd excited state 

mixed 

0 

3 

2 


Spin 

-47r/Lj, 

4 

6 


charge 

Q-nlLy 

3 

6 


mixed 

277 ! Ly 

6 

4 

3rd excited state 

mixed 

4 

6 

— 277 j Ly 


Spin 

— %77lLy 

7 
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the corresponding species of bosons, and Va^ is the ve¬ 
locity matrix. To diagonalize the above Lagrangian, we 
introduce the charge and (pseudo) spin modes (l>c{s) = 

i4>a ± (f>b)lV2. 

We can now obtain the edge Hamiltonian and the cor¬ 
responding momentum operator [22]: 

O'jr 27r 

H=—{v,L^o + VsL^o), P=—{Ll-Ll), (5) 

Ijy l^y 


with 


FIG. 3: (color online). The entanglement spectra versns 
momentum ky-. (a) charge sector ANa + ANb = 0. (b) charge 
sector ANa + ANb = 1. The simulation is carried on an 
infinite cylinder of width Ly = 4 unit cells {Wy = 16 sites), 
a = 7r/2, A = 0. The dashed arrow-lines denote the edge 
modes. 


Edge modes from entanglement spectra: The ex¬ 
istence of symmetry protected gapless edge modes is one 
of the defining properties of the SPT phase. The BIQH 
state has two counter-propagating edge modes, which can 
be identified as a charge mode that carries charge with 
no (pseudo) spin, and a (pseudo) spin mode that carries 
(pseudo) spin with no charge. Thus, as long as the C/(l) 
symmetry (charge conservation) is preserved, backscat- 
tering between the two edge modes is prohibited [19]. 
The BIQH state can be described by an Abelian Chern- 

Simons theory with the AT-matrix 

Thus the behavior of its edge modes is similar to the ones 
of FQH states [1], whose effective Lagrangian is: 

47r 

where a,/3 = A,B and l/2TTdx4>a gives the density of 


rc(s) 

-^0 


, f; 

m—1 


( 6 ) 


Here, Ly is the length of the ID edge; ANa(b) is the 
change in the particle number of a{b) boson relative to 
the ground state; is the set of non-negative inte¬ 

gers describing oscillator modes. These oscillator modes 
exhibit the well-know 1, 1, 2, 3, ... degeneracy pattern 
[!]■ 

As compared to FQH states with only one chiral mode, 
Eq. (11) shows two counter propagating modes and thus 
the BIQH is a non-chiral phase. Here we briefly mention 
how to understand the spectrum of the edge modes of the 
BIQH state. Basically, we expect Vc ^ Vg (when A = 0, 
this is exact due to the symmetry of the Hamiltonian), 
thus Lq+Lq determines which level (ground state, 1st ex¬ 
cited state, ...) of the edge Hamiltonian the mode belongs 
to. For example (when ANa + ANb = 0)) -|- Lg = 0 

gives the ground state and we have Lq + Lq ^ 0 for the 
excited states. Specifically if Lq = 0, the edge mode 
is a charge mode with positive momentum; if Lq = 0, 
it is a (pseudo) spin mode with negative momentum; if 
Lq,Lq yb 0, it is a mixed mode with both charge and 
(pseudo) spin. Once the value of Lq, Lq is fixed, we have 
the freedom to choose different combination of ANa^b) 
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FIG. 4: (color online) Numerical data on the phase transition, 
the cylinder’s width is Wy = 8 sites, (a) Correlation length ^ 
in charge sector 5n — 2 with different number of states kept 
in DMRG simulation, (b) Fidelity (wave-function overlap), 
F{a) = |('0(a — 5a)\il){a + (5a))|. 

and to realize the demanded Lg and Lq. Conse¬ 

quently each mode has certain degeneracies for a given 
energy and momentum. From Eq. (ll)andEq. (12), the 
degeneracies of each mode can be calculated straightfor¬ 
wardly, and is shown in Table I (see the supplementary 
materials for more details). 

Numerically, we can use the entanglement spectra as 
a probe of the edge modes [35]. The numerical results 
from the DMRG simulation are shown in Fig. 3. We 
have plotted two different cases that correspond to the 
C/(l) charge sector AA^o + AA^{, = 0 and AA^-l-AAf, = 1. 
The two counter propagating edge modes, and the mixed 
modes are clearly seen. The counting in each sector from 
our numerics agrees well with the theoretical expectation 
in Table I. If a finite A is turned on so that the (pseudo) 
spin is no longer a good quantum number, the (pseudo) 
spin mode is still robust due to the chiral (anomalous) 
implementation of C/(l) symmetry [10, 11, 18, 19]. This 
is also true in our numerical results, we find that the 
entanglement spectra of finite A remain almost the same 
as A = 0 (Fig. 3). 

Phase transition: We study the phase transition be¬ 
tween the two different BIQH states with varying flux a. 
Interestingly, our numerical results suggest that an in¬ 
finitesimal flux a drives the system into a BIQH states. 
This result is reasonable since any finite flux a breaks 
the time-reversal symmetry, hence favors a BIQH phase. 
Similar physics occurs in Haldane’s honeycomb model 
[38], where an infinitesimal flux stabilizes a Chern insu¬ 
lator. To investigate the properties of the phase tran¬ 
sition, we calculate the correlation lengths and fidelity 
(wave-function overlap) for different fluxes a. The corre¬ 
lation length is calculated using the transfer matrix de¬ 
fined in the infinite DMRG’s algorithm, which determines 
the largest correlation length [33] (see the supplementary 
materials). When the system is deep in the BIQH phase 
[a ^ ±7r/2), the DMRG simulations are fully converges 
and the correlation lengths (Fig. 4(a)) are very small. In 


contrast, near the critical point (a = 0), the correlation 
length diverges as the bond dimension is increased. This 
behavior indicates that the system is near a critical point 
(or phase), where the energy gap is small. In addition we 
measure the fidelity (wave-function overlap) F{a) (Fig. 
4(b)), F{a) = [('(/'(a —(5a)l'0(a-l-(5a))l, as a evolves. The 
fidelity (> 0.9) is very close to 1, thus we can exclude a 
strong first order phase transition (which involves level 
crossing). However, the DMRG simulations are difficult 
near the critical point as the bond dimensions diverges, 
hence a more detailed study (particularly taking account 
of finite size effects) is required to make definite state¬ 
ments about the nature of the transition. 

In general, it is possible that the transition between the 
two opposite BIQH is a direct continuous phase transi¬ 
tion. However, it is important to know whether such 
direct continuous phase transition requires some addi¬ 
tional lattice symmetry ]29], such as the inversion sym¬ 
metry between A, B sublattice of the honeycomb lattice. 
A detailed study of the phase transition and the develop¬ 
ment of a critical theory require more refined numerical 
simulations which are left for a future work. 

Conclusion and Outlook: We have introduced a 
microscopic lattice model that realizes a bosonic in¬ 
teger quantum Hall (BIQH) phase which represents a 
17(1) symmetry protected topological phase. The hard¬ 
core bosonic model is defined on a honeycomb lattice 
with correlated hopping and background flux. Using 
DMRG simulation on an infinite cylinder, we find that 
the ground state shows two characteristic properties of 
BIQH: (i) quantized Hall conductance tr^y = F2 and 
(ii) two counter propagating gapless edge modes. This 
model would be an analog of the Ghern insulator [38] , it 
is of sufficient simplicity to be relevant for optical lattice 
experiments [39, 40] and might furthermore serve as a 
guide to find new physical realization of topological or¬ 
dered states [41]. 

Although we have proved numerically that our model 
has a BIQH ground state, it is very demanding to de¬ 
velop a microscopic theory for it. The correlated hop¬ 
ping term plays an essential role for the emergence of 
the BIQH phase, and it couples the bosons in an inter¬ 
esting way (which naively is a mutual flux attachement 
[19]): if there is a boson at site k, the hopping from site 
i to j is the same as normal hopping; if there is no bo¬ 
son at site k, the boson that hops from site i to j gains 
a negative sign. It is interesting to investigate whether 
the correlated hopping term is also useful to engineer 
other interacting SPT phases. A possible extension is 
to realize a BIQH state with a larger Hall conductance 
(such as axy = 4) by assigning a more complicated back¬ 
ground flux in our model. Another intriguing direction 
is to study the topological phase transition between two 
BIQH phases or one BIQH phase to a topologically triv¬ 
ial phase. For example, by adding some repulsion terms, 
it is possible to drive the system into a topological trivial 
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Mott insulating phase; hence one might be able to obtain 
a deconfined phase transition between the BIQH and a 
topological trivial Mott insulator in our model [28, 29]. 

Acknowledgement .—We thank Yohei Fuji for the col¬ 
laboration on the related project. This work was sup¬ 
ported by the Deutsche Forschungsgemeinschaft (DFG) 
through the collaborative research centre SFB 1143. 

Note added .—After our preprint occurred on the arXiv, 
a related work [42] appeared introducing an alternative 
lattice model to stabilize a bosonic integer quantum Hall 
phase. 
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ALGORITHM 


In our DMRG simulations, we use two ways to wrap honeycomb lattice on a cylinder, as shown in Fig. 5. We 
choose open boundary conditions along the x direction, and periodic boundary conditions along the y direction (with 
identical sites labeled by the same letter). 



FIG. 5: Two cylinder geometries that we use in the numerical simulations. The red thick line is the one-dimensional path we 
choose in the DMRG simulation. 

To calculate the correlation length, we define the transfer matrix T as in Fig. 6 (a). Then the correlation length 
of the charge-0 sector is dehned by the first and second largest eigenvalue Ai ,2 (here Ai is normalized to 1) of the 
transfer matrix T, 

C5n=o = —l/lnA2. (7) 

The correlation in the charge-(5n sector is defined similarly, with the quantum number (charge) shifted by 5n in one 
leg of the transfer matrix, as shown in Fig. 6(b), to calculate the corresponding dominant eigenvalue \sn- Then the 
corresponding correlation length is 

C<5n = — l/lnA^ri,. (8) 



FIG. 6: (a) Definition of transfer matrix and correlation length in charge-0 sector, (b) Definition of correlation length 

in charge-5n sector. 

The correlation length determines the largest correlation length in the infinite cylinder [33], for all possible 
operator-operator correlation function: 

(6I(0)62(r))^e-’-/«, ^<^su. (9) 

Where the operator 6| creates charge 6n, and O 2 annihilates charge 6n. For example bja'-bk creates 1 charge, bjajjbkai 
creates 0 charge. 
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Therefore, instead of calculating various correlation functions, one can simply calculate this single quantity ^sn to 
know the length scale of the largest possible correlation length. 


CONVERGENCE OF DMRG 


To see the whether our DMRG simulations converge, we plot the entropy S and correlation length ^ under different 
system size, in Fig. 7, for the system deep in the BIQH phase {a = 7 r/ 2 ). From the entropy, we find that for 
Wy = 8 ,12, simulations have fully converged as we keep 6000 states. Specihcally, the entropy difference as we keep 
4000 and 6000 states is. 


5(6000) - 5(4000) 
5(6000) 


0 . 0002 , Wy = 8 ; 


5(6000) - 5(4000) 
5(6000) 


0.009, Wy = 12. 


( 10 ) 


For Wy = 16, the simulation is not fully converged, where (5(6000) — 5(4000))/5(6000) « 0.03. But it’s already good 
enough for us to extract many physical properties. 

From the correlation length, we can see that ground state is gapped, since the correlation length is small for all the 
system sizes, and it decreases as the system size increases. 



FIG. 7: (a) Entropy versus different number of states kept. 

Here we show Wy = 8,12,16. 



(b) Correlation length versus different number of states kept. 


THE COUNTING OF THE EDGE MODES IN THE BIQH STATE 

The edge Hamiltonian and the corresponding momentum operator of a BIQH state is [22]: 


Ptt 

H = —iv,L-,+VyL^o), 

Ly 


Ott 

p= 

Ly 


with 


( 11 ) 





( 12 ) 


Here we show how to obtain the degeneracies of each edge mode, according to Eq. 11, 12. For example, the spin 
mode of the 3rd excited state (in AiVa + AiVf, = 0 sector) corresponds to Lg = 0 and Lg = 3. Then we have two 
possibilities, AV^ = AVf, = 0 or AVq — AA;, = ± 2 : 

1. AAq = AA(, = 0, thus ^ mn^ = 3. Correspondingly, 

• nf = 3, and = 0 if m 7 ^ 1 
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• nf = 1 , n| = 1 , and = 0 if m ^ 1 , 2 

• n| = 1, and = 0 if m 7 ^ 3 

This gives 3-fold degeneracy. This is why a single chiral mode (free bosonic oscillator) shows the well known 1, 
1, 2, 3, • • • degeneracy pattern. 

2. ANa — ANt = ±2, thus ^ = 2. Similar as above, the bosonic oscillator contributes a 2 fold degeneracy; 

together with the 2-fold degeneracy from we have a 2 x 2 = 4 fold degeneracy. 

Thus the spin mode of the 3rd excited state has 7 fold degeneracy. Similarly, we obtain the results for the sector 
ANa + ANt, = 0 shown in Table II, and for the sector ANa + ANt, = 1 shown in Table III. The counting rule for 
the edge mode is very complicated (although fundamentally straightforward), and our numerical data agrees with the 
counting rule. 


TABLE II: The energy levels of entanglement Hamiltonian (edge mode): ANa + ANb = 0. 


Levels 

mode 

ky 


Tg, Tg 

ANa, ANb 

{«-m} 


Degeneracy 

Ground state 

- 

0 

To 


ANa 

= ANb = 0 

= 0 

u™ = 0 

1 

1st excited state 

Lg + Lg = 1 

charge 

“l-K ! Ly 

Lg = 

= 1, Lg = 0 

ANa 

= ANb = 0 

51 ; mng, = 1 

flm = 0 

1 

Spin 

— 2'RlLy 

LS = 

= 0, Lq = 1 

ANa- 

ANa 

- ANb = ±2 

= ANb = 0 

ngj = 0 

ngj = 0 

= 0 

51 ; mn%, = 1 

3 


charge 

'^'K/Ly 

Lg = 

= 2, Lg = 0 

ANa 

= ANb = 0 

51 ; mn%, = 2 

= 0 

2 

2nd excited state 

L§ + Lg = 2 

mixed 

0 

Tg 

= Lg = l 

ANa 

ANa- 

= ANb = 0 

- ANb = ±2 

51 ; mng, = 1 
51 ; mn'ia = 1 

51 ; mn%, = 1 
ng^ = 0 

3 

Spin 

! Ly 

LI-- 

= 0, Lg = 2 

ANa 

ANa- 

= ANb = 0 

- ANb = ±2 

ng, = 0 

ngj = 0 

51 ; mng, = 2 
51 ; mng, = 1 

4 


charge 

Q-nlLy 

Tg = 

= 3, Lg = 0 

ANa 

= ANb = 0 

51 ; mng, = 3 

ngj = 0 

3 

3rd excited state 

L§ + Lg = 3 

mixed 

21^ ! Ly 

Tg = 

= 2, Lg = 1 

ANa 

ANa- 

= = 0 

— A7V;, = ±2 

51 ; mn%, = 2 

51 ; mng, = 2 

X; mng, = 1 

ngj = 0 

6 

mixed 

— 2'K/Ly 

Lg = 

= 1, Lg = 2 

ANa 

ANa- 

= AA^6 = 0 

— A7V;, = ±2 

51 ; mng, = 1 
51 ; mng, = 1 

X; mng, = 2 
51 ; mng, = 1 

4 


Spin 

— Q-K/Ly 

Tg = 

= 0, Lg = 3 

ANa 

ANa- 

= ANb = 0 

— A7V;, = zb2 

ngj = 0 

= 0 

X; mnla = 3 
51 ; mng, = 2 

7 
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TABLE III: The energy levels of entanglement Hamiltonian (edge mode): ANa + ANb — 1. 


Levels 

mode 

ky 

LI, LI 

ANa, ANb 


{«m} 

Degeneracy 

Ground state 

- 

0 

L§ = Lg = 1/4 

ANa 

- ANb = ±1 

= 0 

ng, = 0 

2 

1 st excited state 

charge 

2^ / Ly 

Lg = 5/4, Lg = 1/4 

ANa 

- ANb = ±1 

I] mn/a = 1 

ng, = 0 

2 

L§ + LS = 3/2 

Spin 

— 2'K/Ly 

Lg = 1/4, Lg = 5/4 

ANa 

- ANb = ±1 

n™ = 0 

5] mn/a = 1 

2 


charge 

^Tv/Ly 

Lg = 9/4, Lg = 1/4 

ANa 

- ANb = ±1 

5] mn/a = 2 

ng, = 0 

4 

2 nd excited state 

mixed 

0 

Lg = Lg = 5/4 

ANa 

- ANb = ±1 

5 ] mng, = 1 

mn/a = 1 

2 

L§ + Lg = 5/2 

Spin 

— drr/Z/j, 

Lg = 1/4, Lg = 9/4 

ANa 

ANa 

- ANb = ±1 

- ANb = ±3 

= 0 

n/n = 0 

5 ] mng, = 2 
ng, = 0 

6 


charge 

G-KjLy 

Lg = 13/4, Lg = 1/4 

ANa 

- ANb = ±1 

rnn/a = 3 

ng, = 0 

6 


mixed 

2'K/Ly 

Lg = 9/4, Lg = 5/4 

ANa 

- ANb = ±1 

mn/a = 2 

5 ] mng, = 1 

4 

3rd excited state 

-^'0 + -^'0 = 7/2 

mixed 

—2KlLy 

Lg = 5/4, Lg = 9/4 

ANa 

ANa 

- ANb = ±1 

— ANb = i3 

5 ] mn/a = 1 
mn/a = 1 

mng, = 2 

ng, = 0 

6 


Spin 

— GK/Ly 

Lg = 1/4, Lg = 13/4 

ANa 

ANa 

- ANb = ±1 

— ANb = i3 

= 0 

= 0 

5] mng, = 3 
5 ] mng, = 1 
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